Search for: All records

ABSTRACT One-dimensional (1D) stellar evolution models are widely used across various astrophysical fields, however they are still dominated by important uncertainties that deeply affect their predictive power. Among those, the merging of independent convective regions is a poorly understood phenomenon predicted by some 1D models but whose occurrence and impact in real stars remain very uncertain. Being an intrinsically multi-D phenomenon, it is challenging to predict the exact behaviour of shell mergers with 1D models. In this work, we conduct a detailed investigation of a multiple shell merging event in a 20 M$$_\odot$$ star using 3D hydrodynamic simulations. Making use of the active tracers for composition and the nuclear network included in the 3D model, we study the merging not only from a dynamical standpoint but also considering its nucleosynthesis and energy generation. Our simulations confirm the occurrence of the merging also in 3D, but reveal significant differences from the 1D case. Specifically, we identify entrainment and the erosion of stable regions as the main mechanisms that drive the merging, we predict much faster convective velocities compared to the mixing-length theory velocities, and observe multiple burning phases within the same merged shell, with important effects for the chemical composition of the star, which presents a strongly asymmetric (dipolar) distribution. We expect that these differences will have important effects on the final structure of massive stars and thus their final collapse dynamics and possible supernova explosion, subsequently affecting the resulting nucleosynthesis and remnant. 

ABSTRACT The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever-increasing computing power makes now possible to simulate in three-dimensional (3D) part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, and abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same time-scale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur. 

ABSTRACT Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multidimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M⊙ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution; thus, revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions, and compact remnants. 

Abstract We present the state-of-the-art single-zone nuclear reaction networkWinNet, which is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method, along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. With this publication,WinNetwill be publicly available and open source at GitHub and Zenodo. 

ABSTRACT The origin of the proton-rich trans-iron isotopes in the Solar system is still uncertain. Single-degenerate thermonuclear supernovae (SNIa) with n-capture nucleosynthesis seeds assembled in the external layers of the progenitor’s rapidly accreting white dwarf (RAWD) phase may produce these isotopes. We calculate the stellar structure of the accretion phase of five white dwarf (WD) models with initial masses ≥ 0.85 $$\, \mathrm{M}_\odot$$ using the stellar code mesa The near-surface layers of the 1, 1.26, 1.32 and 1.38 $$\, \mathrm{M}_\odot$$ models are most representative of the regions in which the bulk of the p nuclei are produced during SNIa explosions, and for these models we also calculate the neutron-capture nucleosynthesis in the external layers. Contrary to previous RAWD models at lower mass, we find that the H-shell flashes are the main site of n-capture nucleosynthesis. We find high neutron densities up to several 1015 cm−3 in the most massive WDs. Through the recurrence of the H-shell flashes these intermediate neutron densities can be sustained effectively for a long time leading to high-neutron exposures with a strong production up to Pb. Both the neutron density and the neutron exposure increase with increasing the mass of the accreting WD. Finally, the SNIa nucleosynthesis is calculated using the obtained abundances as seeds. We obtain solar to supersolar abundances for p-nuclei with A > 96. Our models show that SNIa are a viable p-process production site. 

Meteoritic analysis demonstrates that radioactive nuclei heavier than iron were present in the early Solar System. Among them, 129I and 247Cm both have a rapid neutron-capture process (r process) origin and decay on the same timescale (≃ 15.6 Myr). We show that the 129I/247Cm abundance ratio in the early Solar System (438±184) is immune to galactic evolution uncertainties and represents the first direct observational constraint for the properties of the last r-process event that polluted the pre-solar nebula. We investigate the physical conditions of this event using nucleosynthesis calculations and demonstrate that moderately neutron-rich ejecta can produce the observed ratio. We conclude that a dominant contribution by exceedingly neutron-rich ejecta is highly disfavoured. 

Asymptotic giant branch stars are responsible for the production of most of the heavy isotopes beyond Sr observed in the solar system. Among them, isotopes shielded from the $r$-process contribution by their stable isobars are defined as $s$-only nuclei. For a long time the abundance of ${}^{204}\mathrm{Pb}$, the heaviest $s$-only isotope, has been a topic of debate because state-of-the-art stellar models appeared to systematically underestimate its solar abundance. Besides the impact of uncertainties from stellar models and galactic chemical evolution simulations, this discrepancy was further obscured by rather divergent theoretical estimates for the neutron capture cross section of its radioactive precursor in the neutron-capture flow, ${}^{204}\mathrm{Tl}$( $t_{1/2}=3.78\mathrm{yr}$), and by the lack of experimental data on this reaction. We present the first ever neutron capture measurement on ${}^{204}\mathrm{Tl}$, conducted at the CERN neutron time-of-flight facility n_TOF, employing a sample of only 9 mg of ${}^{204}\mathrm{Tl}$produced at the Institute Laue Langevin high flux reactor. By complementing our new results with semiempirical calculations we obtained, at the $s$-process temperatures of $kT\approx 8\mathrm{keV}$and $kT\approx 30\mathrm{keV}$, Maxwellian-averaged cross sections (MACS) of 580(168) mb and 260(90) mb, respectively. These figures are about 3% lower and 20% higher than the corresponding values widely used in astrophysical calculations, which were based only on theoretical calculations. By using the new ${}^{204}\mathrm{Tl}$MACS, the uncertainty arising from the ${}^{204}\mathrm{Tl}(\mathrm{n},\gamma )$cross section on the $s$-process abundance of ${}^{204}\mathrm{Pb}$has been reduced from $\sim 30\%$down to $+8\%/-6\%$, and the $s$-process calculations are in agreement with the latest solar system abundance of ${}^{204}\mathrm{Pb}$reported by K. Lodders in 2021. Published by the American Physical Society2024